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Abstract. We examine a model in which a nonequilibrium phase transition from an 
active to an extinct state is observed. The order of this phase transition has been 
shown to be either continuous or first-order, depending on the parameter values and 
the dimension of the system. Using increasingly large clusters, we use the cluster 
approximation method to obtain estimates for the critical points in 1-|-1 dimensions. 
For the continuous phase transitions only, extrapolations of these approximations show 
excellent agreement with simulation results. Further, the approximations suggest that, 
consistent with simulation results, in 1+1 dimensions no first-order phase transitions 
are observed. 



1. Introduction 

The order of phase transitions in 1-1-1 dimensions has long been a topic of debate among 
researches in the field of nonequilibrium phase transitions. While it has been argued 
that first-order phase transitions are impossible in this dimension (see for example [1]), 
most agree that such transitions are feasible. For example, Dickman and Tome sought to 
find the simplest model with short-range interactions that exhibited a first-order phase 
transition in one spatial dimension p|. They examined the pair- and triplet-creation 
models with the reactions 

nA — ^ {n + 1)A and A — > (1) 

with n = 2, 3 respectively. The mean field (MF) of such reactions yields a first-order 
phase transition, yet often, continuous transitions are thought to be observed in (1-1-1)- 
dimensional monte carlo (MC) simulations. They varied the diffusion rate D to see what 
effect this had since they expected that, with larger diffusion, the model would exhibit 
more MF-like behaviour due the better mixing of particles, and hence, a first-order phase 
transition might be observed. They found a first-order phase transition for the triplet- 
creation model only. A continuous phase transition was observed in the pair-creation 
model even when 95% of the attempted moves were diffusive. For the triplet-creation 
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model, hysteresis was observed, indicative of a first-order transition, for sufficiently high 
diffusion rate. Recently, Fiore and de Oliveira [3] used the conservative diffusive contact 
process, in which the number of particles is fixed, to confirm Dickman and Tome's 
original findings (see also [Hl^)- Cardozo and Fontanari [5] also examined the triplet- 
creation model and found that the critical exponents change continuously from DP 
values to compact directed percolation values as the diffusion rate was increased. Due 
to the strong crossover effects, however, they were unable to locate the precise position 
of the point where the transition changed order. 

In this paper, we wish to examine a slightly modified version of the pair-creation 
model in the (l-l-l)-dimensional case through both simulation and, importantly, MF 
techniques. The model has previously been shown, by MF, to exhibit both continuous 
and first-order phase transitions whose lines in parameter space meet at a tricritical 
point [6] (see also [7l[Hll9]|). In (l-l-l)-dimensional MC simulations, however, the model 
exhibits a continuous phase transition across the whole phase space. To examine this 
model further, we employ a technique originally introduced by ben-Avraham and Kohler 
[To] which they called the n-site cluster approximation. In most cases, as n increases, 
the method predicts increasingly accurate behaviour of the model in question. In this 
paper, we collect data for n < 5 and extrapolate our findings as n ^ oo to obtain 
approximations for the order of the phase transition and value of the critical point. 



2. The model 



We have a d-dimensional square lattice of linear length L where each site is either 
occupied by a single particle or is empty. A site is chosen at random. The particle on 
an occupied site dies with probability pd, leaving the site empty. If the particle does 
not die, a nearest neighbour site is randomly chosen. If the neighbouring site is empty 
the particle moves there and produces a new individual at the site that it has just left 
with probability k. If the chosen site is, however, occupied the particle reproduces with 
probability pb producing a new particle on another randomly selected neighbouring site, 
conditional on that site being empty. A time step is defined as the number of lattice 
sites N = L'^ and periodic boundary conditions are used. 

We have the following reactions for a particle A for proliferation and annihilation 
respectively, 

A + A + (p — > 3A, A + (p — >2A and A — (2) 

Assuming the particles are spaced homogeneously, the MF equation for the density 
of active sites pit) is given by 

^ = (1 - Pd) pitf (1 - pit)) + Ml - Pd)p(t) (1 - pit)) - PdPit).i3) 

The first two terms consider the sexual and asexual reproduction reactions respectively 
and the final term death of an individual. Equation has three stationary states: 

Po = 0, (4) 
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Figure 1. Mean field phase diagram showing the continuous phase transition occurring 
for /c > Pb and the first-order phase transition for fc < pb- Equations ([MT]) have been 
re-arranged in the figure to make k the dependent variable. 
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For A; > Pb, P+ ^ continuously as Pd 
transition with critical point 

PdX^>Ph) = ^^- 

For k < Pb, we have a jump in p± from (pb — /c)/2pb to zero, this time at the critical 
point 

(k + pbf 



PdSk <Ph)= . , , 

4pb + (fc +Pb)^ 

Further, for k < p-^,, we have a region 

k . {k + pbf 



< Pd < 



(7) 



(8) 



1 + k ^ ^" - 4p^ + (^k + pk)^ 
where the survival of the population is dependent on the population density. In fact, 
we have extinction for 

p{t) < p-{k <pk,Pd). (9) 

For A; < Pb we therefore have a first-order phase transition. The two phase transition 
lines meet at the point = pb, defining the position of the tricritical point k*. At the 
MF level then, we have a phase diagram as shown in figure [TJ In the region to the left 
of the transition lines, there exists at least one real and positive steady state. In the 
shaded region only, there exist two such steady states with none existing to the right 
of the transition lines. The tricritical point not only marks the intersection of the first- 
order and continuous phase lines, but also the line bordering the region with population 
density dependence. 
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3. The cluster approximation technique 

The MF technique that we employed in the previous section assumed independence 
between individual sites. In reality there exist, of course, correlations between nearby 
sites. Increasing improvements to our original MF equation could then be made by 
considering pairs, triplets, quadruplets and so on, of adjacent sites. Such an approach 
was developed by ben-Avraham and Kohler, called the (n, m)-cluster approximation 
method fiO\ . 

We direct the interested reader to ben-Avraham and Kohler's paper [10] and others 
(for example [HI |12] and references therein) for a detailed explanation of the approach. 
Briefly, the method involves examining clusters of size n. For sites 1, 2, ... ,j, we denote 
Psis2...sj(t) as the probability that the j sites are in state {siS2 ■ ■ -Sj} at time t. The 
method is concerned with how PsiS2...s„{t) changes in time and involves deriving a system 
of master equations for P|s-}(t). Since Si = 0,1 in our case, there are 2" such equations. 
It is, however, easy to see that the number of independent equations is much smaller 
than this (see, for example, [TO]). 

Since we are concerned with clusters of size n, we need a way to approximate 
PsiS2-sn from clusters of size N > n. In their {n, m)-approximation, ben-Avraham and 
Kohler consider adjacent clusters of size n with an overlap of m < n sites. Using the 
Bayesian extension process, the approximation is then. 



llj=0 ^s. 

nN—n p 
j=n~m Sj+l---Sj+m 



P _ '■l-j = + + „ , ^ 

-^siS2..-SN T-(N-n T-> ■ ^ ^ 



For example, the (3, 2) approximation for a cluster of size six is given by 

P P P 

p p S2S3S4 S3S4S5 54^5^6 ni"! 

-'S1S2S3S4S5S6 ~ p p p ' \^^) 

S3S4 S4S5 

Ben-Avraham and Kohler found that the (n, n — l)-approximation yields the most 
accurate results and is termed the n-site approximation for short. 



3.1. The 2-site approximation 

The simple 1-site approximation is just our original MF equations, so we examine 
the 2-site approximation. Introducing the subscripts • and o for occupied and empty 
sites respectively, we have two independent variables, chosen to be the particle density 
p(t) = P,{t) and the pair density c{t) = P„{t). It is then simple to derive the other 
2-site probabilities 

d{t) = P.o(t) = pit) - c{t) = Po.it), (12) 
e{t) =Poo{t) = l-2p{t)+c{t). (13) 

To obtain the master equations, we consider the reactions which change the number 
of occupied sites n, and the number of pairs of sites n,,. The reactions are listed in table 
[H along with their probabilities of occurring, given the particle configuration. From this 
table, the master equations for p and c are then derived in the usual way. 
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Table 1. Reactions for the 2-site approximation where, reactions such as • • o» — ¥ 
• o •• for which An, — An,, = 0, have been ignored. A symmetry factor arising 
from the parity symmetry has been included (right column) rather than writing both 
equations down. 



We note that whereas diffusion of the particles did not feature at all in the 1-site 
approximation, it does appear in the 2-site approximation since n,, can be affected (see 
rows four and five in tabled]). 

Deriving the master equations for p and c, we have 

jg^(l-p,) (^-'^'to-''+'^P' -p,p, (14) 

at p 

dc_ (p - c) [(1 - fc) (p^ - c) + (1 - c) (fcp + pbc)] 

^-(1 Pdj ^^^^ 2pdC, (15j 

where we have used such relations as d + e = 1 — p. We notice that, if we make the 
assumption that all the sites are independent so that c = p^, we return to our original 
MF equation. Solving the p-equation, we have the steady states 

Po = 0, 



Ph 



k±Mp^-kY-Ap^[-^^-k 



P± = — 7 -■ (16) 

\ i-Pd J 

These roots are very similar to the roots found in the original MF approximation except 
for the extra prefactor 

Similar to before then, so long as c > 0, for k < pij,we have a first order phase transition 
at 

(^k -\- Ph)^ 

= Ap,+ik+p,r ^'^^ 

We note that critical point p^^ = k/{l + k) is no longer valid because of the denominator 
in the prefactor. To find the critical points then, we solve these equations numerically 
and plot the results for the 1-site and 2-site approximations in figure [21 



Accuracy of the cluster- approximation method in a nonequilibrium model 



6 



k 



0.2 



0.8 



0.6 



0.4 



0, 








0.1 



0.2 



0.3 



0.4 



0.5 



Pd 



Figure 2. Phase diagram according to the cluster approximation method. The hnes 
show the original MF and the markers the numerically obtained values for the 2- 
site approximation. The region in between the two lines and between the different 
markers show the density dependence region. The inset shows the intersection of the 
two markers for the 2-site approximation showing the tricritical point. The horizontal 
hashed lines in both plots show the values of k* according to both approximations. 



To determine the position of the critical points, for values of k and pd, we 
numerically found all of the real steady states with < c, p < 1. Counting the number 
of such steady states indicated which region of the phase diagram we were in. Using 
an iterative procedure enabled us to locate the boundaries between the regions with 
zero, one and two such steady states as outlined in the phase diagram in figure [H For 
increasing p^? the first-order transition lines were indicated by the number of such steady 
states changing from two to zero and from one to zero for the continuous transition. 
The tricritical point k* is then given by the intersection of these two lines. 

We can further test the stability of the steady states by examining the stability 
matrix 



evaluated at the steady states. By calculating the eigenvalues of A, we have that if both 
eigenvalues are negative, the steady state is stable, otherwise it is unstable. Testing the 
non-zero steady states, we find an unstable state in the density dependent region only. 

As we see in figure [21 for small values of k, the analytical value for the critical 
points from the 1-site and 2-site approximations are identical. We see, further, that the 
position of the tricritical point (intersection of the different markers in the figure) is at 
a lower value of k than the 1-site approximation. Numerically, we found it to be at 
k = 0.2139 - less than half of the original MF prediction. 




(19) 
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Figure 3. a) Numerical results for the critical point for various values of k. The 
line shows the original MF approximation (n — 1 and the markers (from right to left) 
the n = 2, 3, 4 and 5. The circles show the numerical simulation results, c) The 
approximation for for k — \ for the different values of n. The circle shows the 
numerical value with the hashed line showing an extrapolation through the points for 
n = A and n — b. The inset shows the tricritical point k*{n) with, again, a hashed line 
showing an extrapolation. 



3.2. Higher order approximations 

Until now, isolated particles have not been considered appropriately. If we aim to build 
a more accurate approximation, we need to be able to consider such particles. By 
examining the 3-site approximation, for example, we consider such probabilities as Po.o- 
Clearly, the price that we pay is an increase in complexity since the number of variables 
and equations increase rapidly. We derived equations for higher order approximations 
for n < 5. The = 5 case required 13 independent variables with over 1,100 reactions 
having to be considered. 

The results for the approximations and simulation results for the critical points are 
shown in figure [3] a) . We clearly see how increasing the size of the clusters gives more 
accurate approximations for the behaviour of the model when comparing to the MC 
simulation results. In particular, for = 1, we show in the main plot of figure [3] b), the 
approximated values for p^^ {n) against 1/n. An extrapolation of the results a.s n —>■ oo 
shows excellent agreement with the MC value. Unfortunately, at the first-order phase 
transition, when the number of real and positive steady state solutions decrease from 
two to zero, such an extrapolation does not lead to good agreement. This is likely to be 
due to the fact that at the continuous phase transition, the correlation length is infinite 
and therefore considering increasingly large numbers of adjacent sites will lead to more 
accurate approximations. At first-order phase transitions, since the correlation length 
remains finite, for n > we would expect the approximations for to be independent 
of n. 

We can further plot the value of the tricritical point as a function of cluster size n 
and, again, extrapolate. As the inset of figure [3] b) shows, the position of k*{n) decreases 
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with n, appearing to become zero for some finite n. In other words, for sufficiently large 
n, the corresponding ra-site cluster approximation would predict no first-order phase 
transition. This, however, would be surprising since we would expect such behaviour 
to be true only as n ^ oo. Since the tri critical transition is infinitesimally close to a 
first-order transition though, we may well expect a tailing-off of this apparent linear 
behaviour for large n. 

4. Conclusions 

We have seen how, for this model, the analytical cluster-approximation method seems 
to correctly predict a continuous phase transition in 1+1 dimensions across the whole 
phase space as n — > oo. Further, the method predicts, with a high degree of accuracy, 
the critical point for continuous phase transitions only. These findings highlight the 
power of the method for this case and will hopefully lead on to further examination of 
the techniques involved. 
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